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Abstract 

For the James breathers in the K2 — K% — K4 chain and for breathers in the i^-chain, we 
prove numerically that these dynamical objects are not strictly time-periodic. Indeed, for the 
both cases, there exist certain deviations in the vibrational frequencies of the individual particles, 
which certainly exceed the possible numerical errors. We refer to the dynamical objects with such 
properties as quasibreathers. For the iCj-chain, a rigorous investigation of existence and stability 
of the breathers and quasibreathers is presented. In particular, it is proved that they are stable 
up to a certain strength of the intersite part of the potential with respect to its on-site part. We 
conjecture that the main results of this paper are also valid in the general case and, therefore, it 
seems that one must speak about quasibreathers rather than about strictly time-periodic breathers. 
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I. INTRODUCTION 

According to the conventional definition 

HQS 

, discrete breathers are spatially local- 
ized and time-periodic excitations in nonlinear lattices. Because of the space localization, 
different particles vibrate with essentially different amplitudes. On the other hand, it is 
typical for nonlinear systems that frequencies depend on amplitudes of vibrating particles. 
Therefore, it is not obvious how a discrete breather can exist as an exact time-periodic 
dynamical object because, in this case, the particles with considerably different amplitudes 
must vibrate with the same frequency. Surprisingly, we did not find an explicit answer to 
this question in the literature on discrete breathers. 

This paper is devoted to some aspects of the above problem. In Sec. Ell we consider the 
breathes introduced by James in Ref. Q| and arrive at the conclusion that there are some 
deviations of the vibrational frequencies of the individual particles from the average breather 
frequency and these deviations certainly exceed the possible numerical errors. 

On the other hand, the analytical form of the discrete breathers used in is not 

an exact solution to the nonlinear dynamical equations of the K 2 — K 3 — K4 chain and, 
therefore, one can suspect that the above deviations are induced by an inaccuracy of the 
initial conditions for solving the appropriate Cauchy problem. 

To establish results beyond suspicion, we consider discrete breathers in a nonlinear chain 
with a uniform on-site and intersite potential of the forth order (see Sec. In other 

words, we study the K 2 — K3 — K4 chain with K 2 = K3 = and call it " ^-chain" . 
For this case, there exists a localized nonlinear normal mode (NNM) by Rosenberg j^, E| 
which represents an exact discrete breather (DB). Let us note that DBs for such potentials 
were discussed in a number of papers (see, for example, Refs. js,!^), but from a somewhat 
different point of view. Here, we obtain practically exact form of the DB and study its 
stability. It turns out that any infinitesimal vicinity of the exact breather solution consists 
of stable dynamical objects (for the appropriate strength of the intersite potential) which 
are not time-periodic. The strict periodicity occurs only along a certain line in the space 
of possible initial conditions which give rise to NNM. All other initial conditions generate 
the dynamical objects which can be considered as quasibreathers, because they correspond 
to quasiperiodic motion. As a consequence, there are some deviations in the vibrational 
frequencies of the individual breather's particles similar to those for the James breathers. 
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Since, in every physical or computational experiment, we cannot tune exactly onto a line 
in the many- dimensional space of all possible initial shapes of the desired periodic solution, 
it is reasonable to speak only about quasibreathers. It seems that such situation occurs not 
only for the considered case admitting the exact solution, but also for the general case. 

In connection with the above mentioned term "quasibreathers", let us note that the term 
"quasiperiodic breathers" is used in literature for different dynamical objects (see Conclusion 
to this paper). 



II. JAMES BREATHERS 



An approximate analytical form of breathers with small amplitudes for the K 2 — K 3 — K4 
chain was obtained in [4J. Some computational experiments with these breathers were 
presented in j^. The main results of can be outlined as follows. 

Let us consider a nonlinear chain of N = 2N + 1 identical masses (m = 1) which 
are equidistant in the equilibrium state. Interaction only between the nearest neighboring 
particles is assumed. Then dynamical equation for the K 2 — K3 — K4 chain reads 

x n = V'(x n+ i - x n ) - V'(x n - x n -i), n = -N..N, (1) 



V{u) = K 2 u 2 + K 3 u 3 + K A u\ (2) 

Here x n is a displacement of n-th particle from its equilibrium position, V{u) and V'{u) are 
the potential of the interparticle interaction and its derivative, respectively. 

It was proved by James that for any Ub (breather frequency) slightly exceeding the 
maximal phonon frequency u max (in our case, K 2 = \ and, therefore, uj max = 2), i.e. 
uj\ = 4 + /i, /iCl , there exists the following breather solution to Eqs. (jl|2j) : 

vn(t) = l - 1 yM-^M i (3) 

V d cosh \ n^f\x) 

where S = iy( 4 )(0)-(y( 3 )(0)) 2 , 

y n = V (u n ) = 2K 2 u n + 3K 3 ul + 4K 4 ul, u n = x n -x n - 1 . (4) 

Here y n are new variables introduced instead of the old variables x n (actually these new 
variables represent forces acting on the particles of the chain). 
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FIG. 1: Mean square deviation rj(t) for the James breather for long-time interval. All parameters 
are the same as in Table [TJ Time t is given in periods To = ir of the 7r-mode. 

Thus, Eq. Q determines a family of breather solutions. Indeed, there exist a breather 
with amplitude proportional to ^//i for any fixed frequency Uf,. The smaller the deviation of 
the breather frequency u>b from the maximal phonon frequency 2, the less the amplitude of 
the breather. For the case — > 0, the hyperbolic cosine in the denominator of © goes to 
unity for all numbers n and, therefore, the breather localization get worse. Actually, in this 
limit, breather tends to the extended 7r-mode with the infinitesimal amplitude 1181. 

n 

Computational experiments reported in |5J have confirmed the theoretical breather shape 
(J3J). We tried to reproduce some results of that paper, for example, those depicted in Fig. [TJ 
To this end we started with Eq. (JHJ) for t — 0, solved the cubic equations (J3J for obtaining 
the initial conditions for x n (t), and then integrated numerically the differential equations 
of the K% — K3 — K4 chain. 

Using the numerical values of K%, K3, K4, Ub and other parameters from the paper 0, 
we indeed obtained some localized dynamical objects which seemed to be time-periodic at 
first sight. But a closer examination revealed more complexity. 

Bearing in mind the question posed in Introduction, we began to follow the evolution of 
frequencies of the individual particles participating in the breather vibration. Some results 
of this analysis are presented below. 

In Table HI we give the frequencies of nine particles (n = —4.. 4) near the center (n = 0) 



4 



of the breather which were calculated within certain time intervals close to the instants tk 
listed in the first row of the table. Note that all these frequencies u>j(tk) are sufficiently close 
to the breather frequency u h = 2.01, which was used in Eq. (J3J), but their deviations from 
uJb certainly exceed the possible numerical errors. 

Let us comment on the computational procedure. We used the fourth-order Runge-Kutta 
method with time step h of about 0.0003 To, where To = = n. For the times given in 

^max 

Table |U our simulations conserved the total energy of the chain up to 10 -10 . The frequencies 
of the individual particles were obtained by calculating adjacent zeros of the functions Xj(t) 
in certain intervals near fixed instants tk- In turn, these zeros were computed by dichotomy 
and by Newton-Rafson method. 

It is expedient to introduce certain mean values characterizing frequency deviations of 
the individual breather particles. We specify the mean value u(tk) and the mean square 
deviation r](tk) of different ojj{tk) for the breather particles at the moment tk as follows: 

u(t k ) = -^2^j{tk), (5) 



k 
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3 M(M-l) ■ (6) 
Here j = —M..M (M < N) are numbers of the breather particles with significant values of 
Xj(t). The values u>(tk) and r)(tk) are given, respectively, in the two last rows of Table 1. In 
the last column of this table, we give 



/ Et>j(**)-*(**)] 2 

aj V M(M-l) 

which represents the mean square deviation of the frequency for each breather particle after 
averaging upon different moments (for this averaging, we have used all the frequencies u>j(t) 
which were calculated up to t = tk indicated it Table 1). [Note that M in Eqs. (I5I6|) is the 
number of considered breather particles, while M in Eq. (J2J) is the number of u)j(t) which 
were calculated]. 

Some questions arise in connection with the result. 

1. Because f](tk) are rather small, one can suspect that they brought about by certain 
numerical errors. Is it true? 
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TABLE I: Deviations in frequencies of the individual particles for the James breather in K 2 — if 3 — 
K A chain with N = 101 particles. Here K 2 = 0.5, K 3 = -0.1, K 4 = 0.25, uj b = 2.01, /x = 0.0401, 
T = ir,h = 0.001 





ii = 100 T 


t 2 = 600 T 


t 3 = 1000 T 


t A = 1500 T 


a 




2.0101459926 


2.0097251337 


2.0109780929 


2.0085043104 


2.3348301075e-5 


w ■? 

— 


2.0068349042 


2.0093263880 


2.0070499430 


2.0100257555 


2.4686634177e-5 


W 9 
w — z 


2.0112977946 


2.0089800520 


2.0110369763 


2.0080931884 


2.6957224898e-5 


— 1 


2.0065877923 


2.0099366007 


2.0070857334 


2.0101152322 


2.9130307350e-5 




2.0117517055 


2.0086968158 


2.0110115290 


2.0080038280 


2.9739839019e-5 




2.0065877923 


2.0099366007 


2.0070857334 


2.0101152322 


2.9130307350e-5 


U 2 


2.0112977946 


2.0089800520 


2.0110369763 


2.0080931884 


2.6957224898e-5 


^3 


2.0068349042 


2.0093263880 


2.0070499430 


2.0100257555 


2.4686634177e-5 


CJ4 


2.0101459926 


2.0097251337 


2.0109780929 


2.0085043104 


2.3348301075e-5 


U) 


2.0090538525 


2.009403685 


2.00925700227051 


2.0090534223 




77 


7.61439353253e-4 


1.5116218526e-4 


6.9232221345e-4 


3.2690761972e-4 





2. Formula © represents only an approximation to the unknown exact breather solution 

nn 

since it was obtained by neglecting some higher-order terms |5( . To what extent the 
appearance of nonzero rj(tk) reflects the properties of the real breather? 

3. Is there any growth of rj(tk) for large times? It is an important question because 
such growth, if exists, possibly means the onset of stability loss of the exact breather 
solution. 

To shed some light on numerical errors problem, we computed f](tk) and Oj for the ir- 
mode which represents a strictly periodic dynamical regime in any nonlinear chain (see, for 



example. 



111]). Using the same computational procedure as that for obtaining Table HJ 
we get Table ILT1 for the case of the 7r-mode vibrations. From the latter table, one can see that 
deviations rj(tf.) and Gj for the 7r-mode [n| turn out to be zero (up to machine precision). 
Comparing these results with those from Table HI we conclude that deviations in vibrational 
frequencies of the individual particles for the James breather are not numerical errors. 

The function rj(t) is depicted for large time intervals in Fig. ^ This function was found 
by calculating all zeros of each displacement Xj(t) and averaging the obtained frequencies 



TABLE II: Frequencies of the individual particles for the 7r-mode with amplitude A = 0.1. Here 
K 2 = 0.5, K 3 = -0.1, K A = 0.25, T = 7T,h = 0.001 





100 T 


600 T 


1000 T 


1500 T 


a 




2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





U) 3 
w — o 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





9 
— z 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





LA/ — ^ 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 


o 




2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





LJi 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





CJ 2 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 







2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





u;4 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 





(D 


2.00748367227 


2.00748367227 


2.00748367227 


2.00748367227 
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ojj(t) with the aid of Eq. (jfi|). From this figure, it is obvious that rj(t) does not increase in 
magnitude, and demonstrates certain oscillations similar to chaotic. 

The above discussed behavior of the breather particles will be interpreted in the next 
sections with an example of a model which admits an exact breather solution. 

III. EXISTENCE OF BREATHERS IN THE K 4 -CHAIN 

We consider 2N + 1 particles chains with fourth-order potential and periodic boundary 
conditions. The potential includes both on-site and intersite parts. In contrast to the on-site 
terms, the intensity of the intersite terms will be varied. We write this potential in the form 

^iE^+iE^-^-i) 4 . (8) 

n n 

The Newton dynamical equations for such chain can be written as follows: 

x n = -xl + (3[(x n+1 - x n ) 3 - (x n - Xn-i) 3 }, n=-N..N, (9) 

X N+1 (t) = X_ N (t), X_AT_i(t) = X N (t). (10) 
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Here j3 is the parameter characterizing the intensity of the intersite potential. It will be 
shown that stability of the breather solution in the chain @ depends essentially on the 
value of (3 (see below). 

Models similar to @ have been considered in the papers 0,0], but we analyze the chain 
© with different purposes and in a different manner. 

It is well-known that space and time variables can be separated in Eq. © and this was 
done in |8j (for the case without on-site potential). We prefer to treat breather solution 
to Q in terms of the nonlinear normal modes (NNM) introduced by Rosenberg in (| 
Indeed, it was proved that for any uniform potential there exist (localized or/and delocal- 
ized) NNMs which represent strictly periodic motion of all the particles of the considered 
mechanical system. More precisely, equations of motion of a many particle system for the 
dynamical regime corresponding to a fixed NNM reduce to only one differential equation 
for the displacement Xo(t) of an arbitrary chosen particle (this is the so called "leading" or 
"governing" equation), while displacements Xj(t) (j 7^ 0) of all the other particles are pro- 
portional to xo{t) at any instant t. Such dynamical behaviour is remeniscent of the linear 
normal modes whose time dependence is represented by sinusoidal functions because the 
leading equation, in this case, is the equation of the harmonic oscillator. 

It is worth to mention that "bushes of modes" introduced in [3] and investigated in a 



number of other papers 



110 . 
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13. 



14. 



15j | represent a quasiperiodic motion because we have 
m leading differential equations for the m-dimensional bush and, therefore, NNMs should 
be thought of as one- dimensional bushes. 

Below, we consider the procedure for obtaining NNMs. Assuming 

x n (t) = k n x (t), n = -N..N, (11) 

for any time t with constant coefficients k n and substituting this expression into Eq. (JHJ), we 
obtain 

k n x = {-k 3 n + (3[(k n+1 - k n f - (k n - k n ^) 3 }}x 3 , (12) 

where 

kN+i = k_ N , fc_jv_i = k N (13) 

in accordance with the boundary conditions (|1U|). 

The leading equation (it corresponds to n = 0) reads: 

£ = {-1 + 0P1 - 1) 3 - (1 - M 3 ]K (14) 
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because we can assume k = 1. 

Demanding all other equations (fT2"|) to be identical to Eq. (|14jh we obtain the following 
relations between the unknown coefficients k n {n = —N..N, n^O): 

-kl + (3[(k n+ l - k n f - {k n - k^f] = k n {-\ + /3[{h ~ l) 3 - (1 - k^f]}. (15) 

Thus, we arrive at the system of 2N algebraic equations with respect to 2N unknowns 
k n (Eqs. must be taken into account). 

Any solution to Eq. (JT5|l determines a certain form of NNM or its spatial profile. In our 
case, there are some localized and delocalized modes among these solutions. In particular, 
one of the solutions to Eq. (|15|) represents the 7r-mode. Obviously, every localized NNM 
is an exact discrete breather in accordance with its definition as spatially localized and 
time-periodic vibration. The time dependence of the breather is determined by the leading 



the Jacobi elliptic function cn (t, ^ J (see 



equation (j!4j) which can be solved in terms o 
below). Note that this result was obtained in 

In this paper, we will be interested only in the breather solution which is symmetric with 
respect to its center. Therefore, we must assume the following relation to hold: 

k_ n = k n , n = -N..N. (16) 

Taking into account Eq. ()16|) allows us to reduce by a half the number of unknowns in 
Eq. Ijl5j) . Let us write down these equations for the cases N = 1 and iV = 2 in the explicit 
form. 

For N=l, we have the chain with three particles only (iV = 3). In this case, the dynamical 
equations © read: 

'±-\ = -x 3 .! + /3[(x - x_i) 3 - (x-i - Xi) 3 ], 

%o = -%l + ^[(^l - x of - (>o - x-i) 3 ] , ( 17 ) 
Xi = —x 3 + /3[(x-x — Xi) 3 — (xi — Xo) 3 }. 

According to Eqs. (fl"l"|) and ((EH) , the symmetric breather pattern (/c_i — ki = k) reads 

{kx (t), x (t), kx (t)}. (18) 



Substituting this pattern into Eqs. (fTTj) . we o 
the condition of identity of all these equations 



3tain the following algebraic equation from 
3: 



k(l + k)+(3(l-k) 2 (l + 2k) = 0. (19) 



Since the root k = k(/3) of Eq. (fT9"j) is a function of (3, let us choose j3 = 0.3 (below, it will 
be shown that the breather is certainly stable for this value of (3). 
Using MAPLE, we obtain 

k = -0.29344944496399996095. (20) 

Let us note that calculating all the Rosenberg modes with the aid of Eq. (JT5J) we used 
the MAPLE specification Digits=20. Nevertheless, some values that have been calculated 
with such precision will be, for compactness, presented with a smaller number of digits. 

Analogously, for the case N = 2, the chain consists of five particles (N = 5) and we must 
search the symmetric breather pattern as follows 

{k 2 x (t), hxoit), x (t), kix (t) ,k 2 x (t)}. (21) 

Then the algebraic equations for k\ and k 2 read 

h[-l + 2(3(h - l) 3 )] = -kf + (3[(k 2 - hf - (h - l) 3 ], 
k 2 [-l + 2(3(h - l) 3 ] = -k 3 2 - [3(k 2 - 

and we obtain the following roots of these equations for (3 = 0.3: 

fci = -0.29928831163054746768, 

(23) 

k 2 = 0.00359934143244973138. 

Continuing in this manner, we find symmetric breathers for the chains with iV = 
7, 9, 11, 13, etc. particles. Some results of these calculations are presented in Table UTTl 
Being calculated with 20 digits, these results are practically the exact breather solutions for 
the corresponding K4 chains. Moreover, comparing the profiles for the chains with iV = 9 
and N = 15 particles, one can reveal that the further increase of N does not affect the 
spatial profile of the breather solution. Indeed, the considered breathers demonstrate so 
strong localization that the displacements of the particles distant by more than three lattice 
spacings from the breather center are utterly insignificant (they don't exceed 10~ 20 and we 
denote them in the tables by asterisk). Therefore, we conclude that the profile for N = 19 
(and even for N = 5) can be considered as that for the infinite chain (N = 00). 

For the case without the on-site potential we have for N = 3: k = 1, — k± — 
—0.5. The similar results for N = 5,9, 15 are presented in Table HVl From this table, it is 
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TABLE III: Spatial profiles of symmetric breathers in the ^-chain with N = 5,9, 15 particles for 
P = 0.3 





7V" — ^ 

I V O 


AT — Q 


7V" — 1 ^ 

1 V — J. O 


K -7 






* 
* 


h r 

h , 

t\j 4 






* 
* 






£0/101 7/171 /l^9^Q1 7840*=. R 


n R0401 7471 4^9^01 77^1 QQ4o 8 




n oo^qq^/i i /i ^9A/iQ7^i ^ri 9 


00^^QQ^A770R9Q9^^90Q79 

u.uuooyyo4i: / (Uozyzoozuy/ z 


n 00'} c ;QQ'}477089Q9 c i c i90Q79 

u.uuooyyo^ / / uozyzoozuy /z 


7<* , 


9QQ9&R^1 1 fi^0^zL7zlfi7fiR 


9QQ9R&'}1 901 ^0079/170/L 


900988^1 901 ^00794704*30 




1 

J. 


1 


i 




-D 29998831 1 R305474fi7fi8 


-fl 99998831 201 300794704 


-D 29998831 901 3007247041 9 


K 2 


n oo^qq^/i i a ^9A/iQ7^i ^ri 9 


00^QQ^/L77089Q9^90Q79 

u.uuooyyo4i / (Uozyzoozuy/ z 


00^00^/17708909^90079 

u.uuooyyo^ / /uozyzoozuy /z 


k 3 




-0.6040174714525917849e-8 


-0.6040174714525917731994e-8 






* 


* 


h 






* 


k 6 






* 


k 7 






* 



obvious that the results of Ref. ja| don't correspond to the exact solution for the case of the 
infinite chain since the author has used the profile only for TV = 3. Indeed, the breathers 
considered in [^] are determined by the pattern {0, .., 0, 0, — ~, 1, — |, 0, 0, 0}, i.e. all the 
particles outside of the central three-particle domain are assumed to have zero amplitudes 
of oscillation. We cannot be sure why the author of that paper refers to such dynamical 
objects as exact breathers despite they represent only a certain approximation. On the other 
hand, it is evident from Tables ILTT] IIVI that dynamical objects with N > 5 particles can be 
practically considered as exact breathers. 

Now let us continue to study the i\~4-chain with on-site and intersite potential for the 
case j3 = 0.3. 

The time dependence of the breather solution is determined by the leading equation (|14|h 
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TABLE IV: Spatial profiles for symmetric breathers in the .fQ-chain without on-site potential 





N = 5 


N = 9 


N = 15 


k 7 
k r 
k k. 

— 




* 


* 
* 
* 
* 


k s 




-0.17336102462887968846e-5 


-0.1733610246288796884739e-5 


k 2 


0.023048199202046015774 


0.023050209905554654592 


0.02305020990555465459272 


fe i 


-0.52304819920204601577 


-0.52304847629530836653 


-0.5230484762953083665346 


fcn 


1 


1 


1 


fei 


-0.52304819920204601577 


-0.52304847629530836653 


-0.5230484762953083665309 




0.023048199202046015774 


0.023050209905554654592 


0.02305020990555465459218 






-0.17336102462887968846e-5 


-0.1733610246288796884198e-5 


A;4 




* 


* 


fc 5 






* 


fc 6 






* 


£ 7 






* 



For the symmetric breather (k-i = k± = k), it can be written as follows: 

x +p 2 x 3 = 0, (24) 

where 

p 2 = 1 + 2/3(1 - A;) 3 . (25) 

The parameter p = p(N) varies slightly with changing the number N = 2N + 1 of particles 
in the considered chain: 

p 2 (3) = 2.2983734517955912888 
p 2 (5) = 2.3160362289746275590 

(26) 

p 2 (7) = 2.3160362301367966972 
p 2 (9) = 2.3160362301367966972 etc. 

For initial conditions 

x (0) = A), *o(0) = (27) 
12 



the solution to Eq. (J2U (see, for example, |8j) reads 

x (t) = A cn (ut, -Lj , (28) 
where the frequency u is the linear function of the amplitude A$: 

u = pA . (29) 

Here cn(wt, m) is the Jacobi elliptic function with the modulus m equal to Note that 
such value of the modulus is needed to eliminate the linear in Xo(t) term, because, in general 
case, the function cn(r, m) satisfies the equation 2l| : 

cn"(r, m) + [1 — 2m 2 ]cn(r, m) + 2m 2 cn 3 (r, m) = 0. 

Introducing the new time and space variables r, x(r) according to relations 

f=^-, x (t) = Aqx(t), (30) 
PA 



we obtain from Eqs. ITfjl the following Cauchy problem for the function x(t) 

x" + x 3 (t) = 0, x(0) = 1, x\0) = (31) 

with the solution 

x{t) = cn (t, i=J . (32) 

As was the already mentioned, dynamical objects whose existence is derived above demon- 
strate a strong localization and we can describe them using the chain with N = 7 particles 
only. Considering longer chains would not contribute to the accuracy of description. 

The strong localization occurs not only for f3 = 0.3, but also for others /3 (2^ (see 
Tables fV j IVI|1 . As a matter of fact, the localization varies, but this change is completely 
negligible. 

In conclusion, it is worth to emphasize that the spatial profile {/c_ n , &_ n+1 , .., k_ 2 , h-i,ko = 
1, fci, &2> •■■> kn-u k n } provided in Tables ITTT) IIV1 fVl IVII turns out to be universal for the 
breathers with different amplitudes Aq (it does not depend on the amplitude), while the 
breather frequency depends on A linearly (u = pA ). 
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TABLE V: Spatial profiles of symmetric breathers in the .ft^-chain with N = 5,9, 15 particles for 
(3 = 0.5 
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1 V O 
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^3 




-0.83729113342838511470e-7 


-0.83729113342838512537e-7 


/C4 




* 


* 








* 


fc 6 






* 


fc 7 
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IV. STABILITY OF BREATHERS IN THE K 4 -CHAIN 

A. Linearization of the dynamical equations near the breather solution 

To study the stability of a given periodic dynamical regime, in accordance with the stan- 
dard prescription of the linear stability analysis, we must linearize the nonlinear equations 
of motion in the vicinity of this regime (the breather solution, in our case) and investigate 
the resulting linear equations with time-periodic coefficients. 

Let us start our stability analysis with the simplest example, namely, we will consider 
the stability of the breather 

x b (t) = {kx (t), x (t), kx (t)} (33) 
in the three-particle /Q-chain described by dynamical equations (fTTj) . To this end, we 
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TABLE VI: Spatial profiles of symmetric breathers in the ^-chain with N = 5, 9, 15 particles for 
P = 0.6 
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introduce an infinitesimal vector 

8{t) = {6^(t), S (t), S^t)}, 



(34) 



substitute the vector x(t) = x b {t) + 8{t) into Eqs. (|T7|) and linearize these equations with 
respect to Sj (j = —1, 0, 1). 

As the result of this procedure, we obtain the linearized system 



with the symmetric matrix 



A 



S(t) = -3a%(t)A6(t) 



yU + V —jj, ^ 

—fj, 1 + 2/i —fj, 
\ —fj, n + v j 



(35) 



(36) 
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where /jl = (3(1 — k) 2 , v = k 2 . Here the coefficient k is determined by the algebraic equation 
(IT9"j) . while xo(t) is the solution to the leading equation with the initial conditions (|2Tj) 
(for (3 = 0.3 p 2 = p 2 (3) = 2.2983734517955912888). 

It can be easily shown that Eq. (J35|) is valid for an arbitrary value N, but the correspond- 
ing matrix A, in this case, turns out to be more complicated. For example, for the X 4 -chain 
with five particles (N = 5) we have 



' t) + k 



A 



-V 












V 



o 






T) + fi + k\ -fi 

— fJL +1+2^X -yU 

-fi r) + fx + k\ 

-j] 



\ 



-v 

rj + k 2 J 



(37) 



where fi = (3(1 — A4) 2 , r\ = (3(ki — k 2 ) 2 . Here k\ and k 2 are determined by Eqs. [for 
(3 = 0.3, their numerical values are given by Eqs. ()23|) ]. while xo(t) is the solution to Eq. (|2*^ 
with p 2 = p 2 (5) = 2.3160362289746275590. 

The specific structure of the linearized system f!35|) allows us to make an essential step 
in the simplification of our further stability analysis. Indeed, let us pass from the vector 
variable S(t) to a new variable 8(t) whose definition involves a time-independent orthogonal 
matrix S: 

6(t) = S8(t). (38) 

Substituting 6 in such form into Eq. (|3l)j) and multiplying this equation by the matrix S 
from the left (S = S' 1 is the transpose of S), we obtain 24 1 



d = -3x 2 (t)(SAS)d. (39) 

On the other hand, the matrix A is symmetric 25] and, therefore, there exists an orthog- 
onal matrix 5* transforming the matrix A to the fully diagonal form Adi ag '- 



SAS = A, 



diag ■ 



(40) 



If we find such matrix S, the linearized system (J39|) decomposes into N independent differ- 
ential equations 

5j + Sxl^XjSj =0, j = -N..N, (41) 
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where Xj are the eigenvalues of the matrix A. Moreover, solving the eigenproblem Ay=Xy 
for the matrix A, we obtain not only Xj for Eq. (|4ip. but also the explicit form of the matrix 
S from Eq. (J38j) : its columns turn out to be the eigenvectors yj (J = —N..N) of the matrix 
A. 

Each equation of (|41j) represents the linear differential equation with time-periodic co- 
efficient. The most well-known differential equation of this type is the Mathieu equation 

z + [a-2qcos{2t)]z = 0. (42) 

The (a — q) plane for this equation splits into regions of stable and unstable motion [16]. 
If parameters (a,q) fall into a stable region, z(t) that is small at the initial instant t = 
continues to be small for all times t > (the case of Lyapunov stability). In the opposite 
case, if z(0) is a small value (even infinitesimal), z(t) will begin to grow rapidly for t > 
(Lyapunov instability). Actually, we must analyze the stability of the zero solution of Eq. 
(J32J)- In the next subsection, we study the analogous stability properties of the equations 

grj. 



B. Investigation of the basic equation 

Let us consider Eq. (|41jl in more detail. The time-periodic function xo(t), entering this 
equation, is the solution to the Cauchy problem [see Eq. (j2H I27jl ] 

x'q + p 2 x 3 (t) = 0, ar (0) = A , x (0) = 0. 

On the other hand, the change of variables (|3Uj) allows us to eliminate the dependence of 
the Cauchy problem on the breather amplitude Aq (see Eq. (|31|l). Using the same change 
of variables in Eq. (J4"T|) , we get the equation 

4 + ^oW%W = 0, (43) 

where Zj(r) = -^ or ^ e sa ^ e °^ clarity of the further investigation, it is convenient 

to drop the subscript j from this equation and rewrite it in the form 

z" + Ax 2 (t)z(t) = 0, (44) 

U A 3A 

where A = — . 
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FIG. 2: Comparison of the functions cn 2 (r, and |[1 + cos(fir)]. 

Our stability analysis of the breathers in the X 4 -chain is based on Eq. (J44j) and we will 
refer to it as "basic equation". One important conclusion can be deduced immediately from 
this equation, namely, the stability of our breathers does not depend on their amplitudes 
Aq\ The fact is that the amplitude Aq is not contained in Eq. (|44jl neither explicitly, nor 
implicitly (jp is expressed via ki, which, in turn, are determined by an algebraic equation 
independent of A ). 

It is interesting that our basic equation (jUj) can be reduced to the Mathieu equation (J4*2|) 
within a certain approximation. Indeed, taking into account Eq. (|32p. we can write 

..." ^ 1 



x 2 (r) = cn 2 (j-, — j « -(1 + co8(fir)), (45) 

where fi 1.6944. Actually, we approximate the periodic function cn 2 ^r, ^gj by the 
two first terms of its Fourier series. Surprisingly, this fit turns out to be a very good 
approximation, as one can see from Fig. EJ where the function cn 2 (r, -^=) and |(1 + cos(fir)) 
are shown by the solid and dashed lines, respectively. Note that several first terms of the 
Fourier series for cn 2 (t, ^=J read 

cn 2 (t, = 0.4569 + 0.4972 cos(fir) + 0.0429 cos(2fir) + 0.0028 cos(3fir) + .... 
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Introducing a new time variable according to the formula Qt = 2t + tt, we obtain from 
Eq. flSJ): 

2A 

z + —[I - cos(2t)} -z(t) = 0. (46) 



Obviously, this is the Mathieu equation (|42j) with 



2A l A . . 

a=— and q=jp- (47) 



Thus, we obtain a certain relation between two arbitrary parameters of Eq. 

a = 2q. (48) 

Let us now recall the well-known stability diagram for the Mathieu equation (see, for exam- 
ple, [16]). In Fig. El the regions of unstable motion are shaded. 

The boundary lines a,i(q), bi(q) (i = 1,2) of these regions can be approximated by the 
following series: 

2 3 4 

b\ (q) — 1 — o h 

1W y 8 64 1536 ' 

/ x -, , q 2 q 3 q A 
aiiq) = 1 + q -J-64-T536 + -> 

M«) = 4-- + — + ..., 

2W; 12 13824 

In Fig. El we also depict the straight line a = 2q according to Eq. (}4*%|) and four points 

(A,B,C,D) of the intersection of this line with the boundary curves bi(q), a,i{q) (i = 1,2). 
The following values of q correspond to these points of intersection: 

q A = 0.329, q B = 0.890, q c = 1.858, q D = 2.321. 

On the other hand, q — (see Eq. (I47J1 ) and, therefore, we can approximately find the 
values of the parameter A corresponding to the boundaries of stable and unstable motion for 
the basic equation (jUj) . In this way, we obtain the regions represented in the first column 
of Table EIH 

Thus, there exist regions of stable and unstable motion for our basic equation (J44j) in 
accordance with the numeric values of the parameter A entering this equation. 

The above approach based on the Mathieu equation analysis not only sheds light on the 
main properties of Eq. (jUJ), but allows us to arrive at some approximate results presented in 

19 



FIG. 3: The regions of stable and unstable motion for the Mathieu equation. 

Table IVTT1 On the other hand, we can obtain analogous results with high precision analyzing 
the basic equation (jUj) with the aid of the Floquet method. 

To this end, we construct the 2x2 monodromy matrix by integrating Eq. twice (with 
the initial conditions z(0) = 1, z'(0) = and z(0) = 0, z'(0) = 1) over one period of the 
function cn fr, J and calculate its multiplicators. If the absolute value of a multiplicator 
exceeds unity by more then 10~ 5 , we identify the case of unstable motion. 

The results obtained by the Floquet method prove to be surprising! Indeed, all boundary 
values of A are integer numbers, at least up to 1CT 5 : 

A A = 1, A B = 3, A c = 6, A D = 10. 

We suspect that this property of Eq. (jUJ) can be proved exactly, but we could not find a 
proof of this conjecture |2fi| |. 

Thus, for Eq. (jUJ, we obtain the regions of stable and unstable motion presented in the 
second column of Table I VI II 
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TABLE VII: Regions of stable and unstable motion for Eq. Q44[) 



Analysis based, on the Mathieu equation 


Exact Flocjiiet analysis 


Stable or unstable motion 


< A < 0.945 


< A < 1 


stable 


0.945 < A < 2.554 


1 < A < 3 


unstable 


2.554 < A < 5.335 


3 < A < 6 


stable 


5.335 < A < 6.665 


6 < A < 10 


unstable 



C. Quasibreathers in the A^-chains 

The above results for the basic equation f)44j) allow us to make a final step in our breather 

stability analysis. Indeed, we have reduced this analysis to the problem of stability of the 

zero solutions to N individual equations (|43|) . A given breather will be stable, if all the 
3A 

values A 7 - = — ^- (j = —N..N) from Eq. (1431) fall into certain regions of stability motion of 
p 2 

the basic equation (JUJ). 

With high precision, we have computed the eigenvalues Aj of the matrix A [see, for 

3A 

example, Eqs. (|35H37|) ] and the values Aj = — ^ entering the basic equation (J4"4"j) for the 
i^4-chain with different number of particles (N) and different values of the parameter (3 
determining the strength of the intersite potential. It has been revealed that all Aj, with 
exception of Ai, depend considerably on (3 and slightly on N . 

On the other hand, A x = 3 for all N and j3 and this constant coincides, at least up to 
10~ 10 , with the boundary between the first region of unstable and the second region of stable 
motion of the basic equation (J44j) [see Table EH]- Because of the importance of the equality 
Ai = 3, we prove it analytically for N = 3 in Appendix 1. 

According to the Floquet theory, this means that a strictly time-periodic solution cor- 
responds to Ai = 3. Moreover, the eigenvector V\ of the matrix A, corresponding to Ai, 
remarkably coincides with high precision (see also Appendix 1 for the analytical proof of this 
fact for the case N = 3) with the spatial profile of the considered breather for all N and f3. 
Therefore, the infinitesimal perturbation along the vector V\ does not relate to the stability 
of the breather: it leads only to the infinitesimal increasing of the breather's amplitude. 
Thus, studying the breather stability, we must consider only Aj with j / 1. 

In Table rVTH[ for (3 = 0.3, we present the eigenvalues Aj of the matrix A and the 
corresponding values Aj for the K 4 -chains with N = 3, 5, 7, 9 particles. It may be concluded, 
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TABLE VIII: Eigenvalues of the matrix A of the linearized dynamical system for the .fQ-chain 
with different number of particles N {(3 = 0.3) 





N 


= 3 






N 


= 5 








= 2.298373452 




Ai = 3 




Ai = 2 316036234 




Ai = 3 




A 2 


= 0.5880160166 


A 2 


= 0.7675201995 




A 2 = 0.6248090298 


A 2 


= 0.8093254611 




M 


= 0.2934494449 


A 3 


= 0.3830310231 




A 3 = 0.3226215801 


A3 




= 0.4178970630 














A 4 = 0.02626701223 


A 4 


= 0.03402409500 














A 5 = 0.02530830363 


A 5 


= 0.03278226390 






N 


= 7 






N 


= 9 






Ai 


= 2.316036234 




Ai = 3 




Ai = 2.316036234 




Ai = 3 




Ao 


= 0.6248090398 


Ao 


= 0.8093254740 




A 2 = 0.6248090398 


Ao 


= 0.8093254740 




A3 


= 0.3226216094 


As 


= 0.4178971017 




A 3 = 0.3226216094 


A, 


= 0.4178971017 






= 0.02627089232 


Aa 


= 0.03402912093 




A 4 = 0.02627089232 


Aa 


= 0.03402912093 




A 5 


= 0.02531216336 


A 5 


= 0.03278726346 




A 5 = 0.02531216336 


A 5 


= 0.03278726346 




A 6 = 


0.3886030305e - 5 


A 6 = 


0.5033639268e - 


5 


A 6 = 0.3886030305e - 5 


A 6 = 


0.5033639268e - 


5 


A 7 = 


0.3886011392e - 5 


A 7 = 


0.5033614770e - 


5 


A 7 = 0.3886011392e - 5 


A 7 = 


0.5033614770e - 


5 












A 8 = 0.1094511318e - 16 


A 8 = 


0.1417738596e - 


16 












A 9 = 0.1094511318e - 16 


A 9 = 


0.1417738596e - 
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from this table, that for N > 7 only A 2 , A3, A 4 and A 5 are of considerable magnitude and 
correct values of these Aj, at least up to 10~ 5 , can be found already from the case N = 5. 
From Table IVIH we also see that all Aj (j > 1) fall into the first region of stability 
(0 < Aj < 1) and, therefore, the considered breathers are stable for (3 = 0.3. 

In Fig.Hl we present the solution to Eq. (J4~4"|) for A 2 , A 3 , A 4 , A 5 corresponding to a certain 
initial amplitude (the value of this amplitude is inessential because Eq. (}4~4"|) is linear). They 
are not periodic, but stationary solutions in the sense that their amplitudes don't increase 
infinitely in time, as it occurs for the case where Aj fall into an unstable region. Thus, if we 
are in a close (even infinitesimal) vicinity of a given breather, i.e. if all Sj(t) [see Eq. (|41jl] 
are small values at the initial instant t = 0, they continue to be small for all later times 
t > 0. Then according to Eq. (JHEJ), the smallness of the Chebyshev norm ||5j(t)|| c implies 
the smallness of ||<5j(£)|| c . 
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FIG. 4: Solutions to the basic equation (|44[) corresponding to different values of A. 

On the other hand, it follows from Eq. (J3%|) that the solution to the original nonlinear 
equation Q is not periodic in any small vicinity of the exact breather! Indeed, because of the 
relation S(t) = SS(t), each 5j(t) is a certain linear combination of all the 5j(t) (j = —N..N), 
but individual Sj(t) are, in general, quasiperiodic functions. [Even if certain 8j(t) would 
be periodic, their periods are independent of each other and, therefore, the total solution 
^j'(^) U = —N..N) will not be periodic in any case.]. In other words, we arrive at the 
conclusion that arbitrary small vicinity of the exact breather solution consist of quasiperiodic 
solutions which can be naturally called quasibreathers. 

Moreover, in the case of the K 4 -chains with f3 = 0.3, these quasibreathers turn out to be 
stable dynamical objects. Indeed, in a sufficiently small vicinity of the exact breather, the 
quasibreathers are described by the vector x(t) = x b (t) +S(t) [see Eqs. (|3*3*jl and (|3~3jl ] where 
dynamics of S(t) is determined by the linearized system ()35j1 . Above we have demonstrated 
that all || 5j(t) || c are certainly small for any t > 0. Therefore, one can conclude that not 
only the considered breather solution, but also the quasibreather solutions, which are close 
to it, must be stable in the Chebyshev norm. 

Actually, in the .fQ-space, there exist a certain one-dimensional family of the exact 
breathers with different amplitudes Aq = Xq(0) and with the same spatial profile (see, 
Tables ITTT] IIV| fV j IVIjl . A straight line corresponds to this family in the A^-dimensional space 
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of all the conceivable initial conditions Xi(0) (i = —N..N). It is practically impossible to 
tune exactly onto this specific line in the many- dimensional space. 

On the other hand, in any vicinity of this line, we have to deal with the quasibreathers: 
the different particles possess different frequencies and, moreover, these frequencies evolve 
in time. Such a behavior of the individual particles in a quasibreather vibration in the 
i^4-chain can be illustrated by the method used in Sec. 2 for studying the James breathers. 

We investigate the stability of the quasibreather solutions by direct numerical integration 
of the differential equations Q of the considered chain over large time intervals. To this 
end, we choose a certain initial deviation 

6(0) = e ■ {5_ N (0), S o (0), MO)} (49) 

from the exact spatial profile Xb(0) of a given breather, where e is a small parameter, while 
5 n (0) [n = —N..N] are random numbers whose absolute values don't exceed unity. Then 
we solve Eq. (JOJ) with initial condition £c(0) = Xb(0) + 5(0), x(0) = and examine the 
numerical solution x(t) = {x_jv(t), ...,xo(t), ...,cc/v(£)} after a long time. We can scan any 
vicinity of the exact breather by varying e and the random sequence {S n (0) \ n = —N..N} 
from Eq. I n Table IIX1 we present the results of such a calculation for the i-Q-chain with 
N = 7 and (3 = 0.3. We have used the fourth-order Runge-Kutta method with time step 
h = 0.001 and integrated Eq. © up to t = 1500T , where T is the period of the 7r-mode 
(Tq = 7r). The frequencies (uj) of only five breather particles (j = —2. .2) have been taken 
into account, because the vibrational amplitudes of the particles with j = ±3 are very small 
(they are of the order of 10~ 8 ). 

In Table Hxj we also present the mean square deviations r\ = rj(tk) for ujj determined by 
Eq. (JHJ), which have been computed for tk ~ 1500T , and the maximal values of r)(tk) for 
the interval < tk < 1500To. The fact is that f](tk) varies on the considered time interval 
(see, for example, Figs. EJ [HJ) and, therefore, i] max is a more relevant characteristics of the 
deviations in frequencies of the individual particles. 

From Table ITXl one can see the specific quasibreather phenomenon, namely, the deviations 
in the frequencies u>j of the individual particles increases with increasing the parameter e, 
characterizing the deviation of the quasibreather shape from that of the exact breather in 
the Chebishev norm. It is important to emphasize that despite different particles possess 
slightly different frequencies u>j, the quasibreathers are stable dynamical objects. Indeed, we 
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TABLE IX: Derviations in frequencies of different breather (quasibreather) particles 





£ = 


e = 10 -7 


e = 10 -6 


e = 10~ 5 


e = 10~ 4 


£ = 10~ 3 




1 989333684 


1 989394677 


1 989996737 


1 988414794 


1 959673857 




(l) 1 


1 289333684 


1 289333939 


1 28933311 


1 289297098 


1 289077907 


1 287039954 


UJn 


1.289333684 


1.28933432 


1.289333134 


1.289332824 


1.2893981453 


1.2933772905 


U\ 


1.289333684 


1.289333847 


1.289332133 


1.289281901 


1.2890993617 


1.2887949577 


LL>2 


1.289333684 


1.289320432 


1.289203072 


1.288248186 


1.2767350998 




V 





2.8923591607e-6 


2.911839882e-5 


2.397956354e-4 


7.103135734e-3 


1.8891448851e-3 




1.261367077e-10 


1.0051826702e-6 


1.0272244997e-5 


1.002013005e-4 


1.000467033e-2 


2.443460953e-l 
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FIG. 5: Mean square deviations rj = ry(tfc) in frequencies of the breather (quasibreather) particles 
for £ « 10" 9 . Time i is given in periods To. 
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FIG. 6: Mean square deviations r] = rjit^) in frequencies of the breather (quasibreather) particles 
for £ « 10" 6 . Time t is given in periods Tq. 
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did not observe any decay of these objects for e < 1CT 2 up to t = 10 6 T (r} max , even for such 
time interval, does not practically differ from those presented in Table ITXj) . 

However, for e > 1CT 2 , we have observed the decay of the quasibreathers which manifest 
itself in appearance of appreciable vibrational amplitudes of those particles whose amplitudes 
were practically equal to zero in the exact breather solution. 

In conclusion, let us return to Figs. El and HI where we depict rj(tk) as a function of 
the subsequent instants % for which the frequencies Uj were calculated. Sometimes, the 
function ijitf.) demonstrate regular oscillations, sometimes, practically chaotic oscillations. 
Such a behavior can be understood, if one takes into account that the displacement of each 
particle is a superposition of different quasiperiodic functions, as it was shown above in the 
present section. 

D. Stability of breathers with respect to strength of intersite potential 

We consider the AVchain determined by the potential energy (jHJ) inducing the Newton 

dU 

dynamical equations (x n = — — — ) in the form Q. Let us discuss the stability of the 
breathers and quasibreathers in this chain with respect to the parameter (3 which deter- 
mines the strength of the intersite part of the potential energy relative to its on-site part. 
Eigenvalues Xj of the matrix A from Eq. (J3H|) and, therefore, the corresponding values Aj of 
the parameter A entering Eq. (jHJ), depend on (3: Xj = Xj(f3), Aj = Aj{(3). 

In Fig. [7J we present Aj(0) with j = 2,3,4,5 for the A^-chain with A = 3, 5, 7, 9 
particles. These Aj((3) are of the more significant values, as it follows from Table IViHl All 
Aj(f3) which are not depicted in Fig. [7| are small positive numbers. Then one can conclude 
that for (3 G [0,0.554] all Aj((3) remain in the first stability region (0 < A < 1) of the 
basic equation (jHJ). As it has been already discussed this demonstrates the stability of 
the considered breathers (and, therefore, quasibreathers which are close to them) when (3 
increases from zero up to (3 C = 0.554. (Note that our breathers are unstable for (3 < 0). 

On the other hand, we find from Fig. that A 2 intersects the upper boundary (A = 1) of 
the first stability region for N > 3 when f3 exceeds 0.554. This implies the loss of stability 
of the breathers in the considered AVchain. 

Thus, the intersite part of the potential with (3 > must not be too large with respect to 
its on-site part ((3 < 0.554) for breathers (quasibreathers) to be stable. 
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FIG. 7: The functions Aj(/3) for different eigenvalues of the matrix A for the i^-chain with 
N = 3, 5, 7, 9 particles. Note that the plots of A 4 (/3) and A 5 (/3) practically coincide with each 
other (see Table EHH for N = 9). 

Finally, let us point out an one- dimensional subspace of the space of all possible dis- 
placements which becomes unstable when (3 intersects the critical value f3 c = 0.554. This 
subspace is determined by the eigenvector V corresponding to A 2 : 

V = [0.05266431305, -0.7051428715, 0, 0.7051428725, -0.05266431306]. 
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V. CONCLUSION 



The main conclusion of the present paper is that the conventional view on the discrete 
breathers as strictly time periodic and spatially localized dynamical objects must be revised 
in a certain sense. Indeed, it has been shown here that for the James breathers 0, 0] 
in the K2 — K3 — K4 chain as well as for the breathers in the /Q-chain with on-site and 
intersite potentials we actually deal with the dynamical objects representing quasibreathers: 
there are certain deviations in the vibrational frequencies of the individual particles. These 
deviations can be characterized by the mean square deviation which certainly exceeds the 
possible numerical errors. Moreover, for the case of the .fQ-chain, we have performed a 
rigorous investigation of the existence and stability of such quasibreathers. For the K4- 
chain, the exact breathers exist only along a certain line in the many-dimensional space of 
all the possible initial conditions, and it is actually impossible to tune precisely onto this 
line in any physical and computational experiments. 

In some of our numerical experiments with the ZQ-chain, the deviations in frequencies 
of the vibrating particles from the average quasibreather frequency Uj, attained 1%, but 
possibly these deviations can considerably exceed this value for more realistic models. 

The deviations in frequencies of the individual particles (and in frequency of the given 
particle over time) result in some change of the breather Fourier spectrum, namely, instead of 
the ideal lines at the breather frequency Uf, and its multiples there appear certain (possibly, 
narrow) packets of the Fourier components near the ideal breather lines and near zero 
frequency. This effect is difficult to reveal with the aid of the numerical Fourier analysis and 
we prefer to study the deviations in vibrational frequencies of the individual particles in the 
straightforward way. 

It is essential, that the above described deviations in frequencies of the individual vibrat- 
ing particles, in general, don't mean an onset of the breather decay. We have demonstrated 
this fact with the example of stable quasibreathers in the AVchain. For this case, we suc- 
ceeded in proving that these quasibreathers turn out to be stable up to a certain strength of 
the intersite potential with respect to the on-site potential. 

We conjecture that the results obtained in the present paper for two particular cases (the 
James breathers and the breathers in the i^ 4 -chain) are also valid for the general case and, 
therefore, one must speak about quasibreathers rather than about strictly time-periodic 
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breathers. 



Finally, let us note that the term "quasiperiodic breathers" is used in literature for 
dynamical objects different from the quasibreathers considered in the present paper (for 
example, see [itJ and the references therein). Indeed, the former possess several basis 
frequencies (and their integer linear combinations) in the Fourier expansion with substantial 
amplitudes, while the latter possess only one basis frequency u;& (and its multiples), as well 
as many small components with frequencies different from ncob- The quasiperiodic breathers 
exist only in rather specific cases, while quasibreathers seam to be the generic dynamical 
objects. 
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APPENDIX A: 

In Sec. II V C[ using straightforward numerical calculations, we have demonstrated that the 
largest eigenvalue Ai of the matrix A [see Eqs. (jHMHZj) ] corresponds to the boundary between 
the first unstable and the second stable regions for the basic equation ([44)1 . Moreover, the 
eigenvector V\i associated with Ai, determines the direction of the infinitesimal perturbation 
along the considered breather. 

Below we prove these propositions analytically for the JQ-chain with A^ = 3 particles. 
Let us consider Eq. ()35|) 



with matrix A determined by Eq. ([36)1 . The parameter k entering the matrix A determines 
the spatial profile 



6 



3x 2 (t) ■ Ad 



(AI) 



{k,l,k} 



(A2) 



of the breather 



{kx (t),x (t), kx (t)}, 



(A3) 
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while time dependence of the breather particles can be obtained from equations (|24I25|) : 

x +p 2 x 3 = 0, (A4) 

p 2 = 1 + 2/3(1 - A;) 3 . (A5) 

On the other hand, parameter k = k(f3) is a function of the intersite potential stregth f3, 
as it follows from Eqs. (fTTjj) : 

A;(1 + A;)+/3(1-A;) 2 (1 + 2A;) = 0. (A6) 

The eigenvalues Xj (j = 1, 2, 3) of the matrix A can be calculated from the characteristic 
equation of this matrix in the following form: 

MM = \[0- + k2 ) + 3/9(1 - kf] ± i(l - k)y/(l + k) 2 + 2/3(1 - W) + 9/9 2 (l - k)\ (A7) 



Hf3) = k 2 + I3{k-lf. (A8) 

In principle, one can obtain these eigenvalues as explicit functions of the intersite potential 
strength (3 from Eq. (|A6|) and substitute it into Eqs. (jA7|) and ()A8|) . This way is too 
cumbersome and we prefer to use the following trick. Let us express via k from Eq. ()A6j) 
and substitute (3 = (3{k) into Eqs. (|A7[ IA8J) . Then, the square root entering Eq. ()A7j) can 
be explicitly extracted and written in the form 

(1 + A;)(1 + 2A; 2 ) 



1 + 2k 

As a consequence of this extraction, we obtain 



(A9) 



A 3 = -A:. (All) 

The same substitution of /3 = /3(fc) into Eqs. (jA8|) and ()A5|) permits us to write A2 and p 2 as 
follows: 
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2 = l + 2fc 

v 1 + 2k ' v ; 

Comparing Eqs. (|A10|) and (|A13|) . we obtain 

\m=p 2 m. (aw) 

This is an important result since our basic equation (}44j) reads 

z'' + A^(t)x(t) = 0, (A15) 

with 

A = ^r- (A16) 

p 2 {(3) 

3Ai 

Then, we can conclude that Ai = = 3 and, therefore, Ai is a constant with respect 

P 2 {(3) 

to the intersite potential strength (3\ Moreover, as already has been discussed in Sec II V CI 
this value turns out to be the exact boundary between the first region of unstable motion 
(1 < A < 3) and the second region of stable motion (3 < A < 6) of the basic equation (|A15|) . 

Let us recall that numerical calculations have convinced us with high precision that 
Ai = 3 not only for the case = 3, but for any other number of the particles in the K4- 
chain. Moreover, it can be proved that Ai = m for the uniform potential of the arbitrary 
order m (in the present paper, we consider only the case m — 3). 

The eigenvectors Vj (j = 1, 2, 3) of the matrix A corresponding to Xj from Eqs. (|A10HA12|) 
can be written after the substitution (3 from Eq. (jA6|) as follows: 



V = % 1, k],V 2 = [-1, 0, 1], V 3 = [1, -2k, 1], (A17) 

where k = k(/3). We see that V\ = [k, 1, k] represents the vector which coincides exactly 
with the spacial profile (jA2|) of the considered breather. On the other hand, Vj (j = 
1,2,3) are the eigenvectors of the matrix A [see Eq. (jAljl ] of the dynamical equations of 
the iQ-chain linearized near the breather ()A3jl . Therefore, V± = [k,l,k], corresponding to 
A = 3, determines the direction along our breather in the three-dimensional space of all the 
displacements of the particles. In Section HVl using numeric calculations with high precision, 
we have already arrived at the conclusion that this result turns out to be correct not only 
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for 



the case N — 3, but also for the ^-chain with arbitrary number of particles. 
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